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ABSTRACT 

We have investigated the ram-pressure stripping of disk galaxies in clusters at vari- 
ous redshifts and in cluster progenitors; the clusters grow up on a hierarchical clustering 
scenario. We consider a radially infalling galaxy whose initial position and velocity are 
given by a spherical collapse model of structure formation. Moreover, since observa- 
tions show that the intracluster medium (ICM) of nearby clusters is non-gravitationally 
heated, we study the effect of the non-gravitational heating on the ram-pressure strip- 
ping. For a given redshift, we find that ram-pressure stripping has more influence on 
galaxies in more massive clusters. On the other hand, for a given mass, it has more 
influence on galaxies in the clusters at higher redshifts. In particular, we predict that in 
rich clusters at z > 1, most of the galaxies are affected by the ram-pressure stripping. 
While the non-gravitational heating significantly reduces the influence of ram-pressure 
stripping on galaxies in clusters with mass smaller than 1 — 5 x 10 13 M , it does not affect 
the influence in richer clusters. If the ICM is heated non-gravitationally at z >> I, ram- 
pressure stripping does not occur at z ~ 1 — 2 in the progenitors of clusters observed 
at < 2 < 0.5, because the heat makes the ICM fraction of the cluster progenitors 
small. On the other hand, if the ICM is heated non-gravitationally at z ~ for the first 
time, the ram-pressure stripping occurs even at z ~ 3. We compare the results with the 
observations of galaxies in clusters at various redshifts. 

Subject headings: galaxies: clusters: general — galaxies: interactions — galaxies: ISM 
— galaxies: kinematics and dynamics 



1. Introduction 

By comparing observations of galaxies in clusters at z > 0.2 with those in clusters at z ~ 0, we 
infer that some environmental effects in clusters have influences on the evolution of the galaxies. 
Butcher and Oemler (1978, 1984) found that clusters at z > 0.2 have a high fraction of blue galaxies 
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in comparison with nearby clusters, and subsequent works have confirmed this trend (e.g. Lubin 
1996; Rakos and Schombert 1995; Ellingson et al. 2000). Recent observations with the Hubble 
Space Telescope (HST) revealed details of the blue galaxies. Dressier et al. (1994) and Couch et 
al. (1998) found that most of the blue galaxies are normal spirals with active star formation. On 
the other hand, observations have shown that the fraction of SO galaxies decreases rapidly with 
redshift in contrast with the normal spirals (Dressier et al. 1997; Smail et al. 1998; Couch et al. 
1998). These suggest that the blue normal spirals observed in high redshift clusters evolve into 
the non-blue SO galaxies observed in nearby clusters. In fact, observations show that in distant 
clusters there are galaxies in post-starburst phase. (van Dokkum et al. 1998; Poggianti et al. 1999). 
These galaxies may be the ones for which star formation activity is dying down. 

Several mechanisms are proposed that can lead to the color and morphological transformations 
between galaxy classes in clusters, such as galaxy mergers (Bekki 1998), tides by the cluster 
potential (Byrd and Valtonen 1990; Fujita 1998), and tidal interactions between galaxies (Moore 
et al. 1996). One of the strongest candidates is the ram-pressure stripping proposed by Gunn 
and Gott (1972). If a cluster galaxy moves through hot intracluster medium (ICM), it feels ram- 
pressure from the ICM, which is given by picM^rei' wnere Picm is the ICM density and v m \ is 
the relative velocity between the galaxy and the ICM. If the ram-pressure becomes large enough, 
the interstellar medium (ISM) of the galaxy cannot be hold by the gravity of the galaxy and 
is swept away. Numerical simulations demonstrate that in a cluster environment, the stripping 
is likely to occur (Takeda, Nulsen and Fabian 1984; Gaetz, Salpeter and Shaviv 1987; Portnoy, 
Pistinner and Shaviv 1993; Abadi, Moore and Bower 1999; Mori and Burkert 2000). In particular, 
high-resolution three dimensional numerical simulations show that the ram-pressure stripping is so 
effective that it removes 100% of the atomic hydrogen content of luminous galaxies within 10 8 yr 
(Quilis, Moore and Bower 2000). On the other hand, Fujita and Nagashima (1999) investigated 
the influence of ram-pressure on the star formation activity of a disk galaxy. They found that just 
before the atomic hydrogen content is stripped, the star formation rate increases at most a factor 
of 2, but rapidly decreases on a timescale of 10 8 yr after the stripping. After the star formation 
activity, which mainly occurred in the disk, ceases, the galaxy looks like a SO galaxy in both color 
and morphology (Fujita and Nagashima 1999). HI deficient galaxies and galaxies with no strong 
emission-lines seen in cluster cores support the theoretical predictions (e.g. Giovanelli & Haynes 
1985; Dale & Uson 2000). Although ram-pressure stripping alone does not explain the detailed 
morphological features of SO galaxies, such as their large bulge to disk ratios or their conspicuous 
thick disks (Burstein 1979a, b), it may be a principal mechanism of the transformation of spirals 
with active star formation into SO galaxies with inactive star formation. 

However, most of the previous studies dealt with the ram-pressure stripping of a model galaxy 
in a given cluster with arbitrary initial conditions. Moreover, they did not take the evolution 
of cluster structure into account; as will be seen in §2, the structure of high-redshift clusters is 
different from that of nearby clusters even for the same mass. Since it affects the ICM density, 
the velocity of the galaxies, and the efficiency of ram-pressure stripping, it must be considered 
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when we compare the theoretical models with observations of high-redshift clusters. In this paper, 
we investigate ram-pressure stripping in clusters at various redshifts, which grow according to a 
hierarchical clustering scenario; the initial position and velocity of galaxies are given by a spherical 
collapse model of cluster formation. Moreover, since a cluster breaks into smaller progenitors as 
time goes backwards, galaxies in the cluster might have been affected by ram-pressure stripping 
when they were in the progenitors before the present-day cluster formed. Thus, we also consider 
the ram-pressure stripping of galaxies in these progenitors. 

Since ram-pressure is proportional to the density of ICM, the ICM distribution of a cluster 
may be related to the evolution of the cluster galaxies feeling ram-pressure. X-ray observations of 
nearby clusters show that their ICM distributions are generally different from their dark matter 
distributions. In particular, for low temperature clusters, the distributions of ICM are significantly 
flatter than those of dark matter and the ICM fraction in their central regions is small (Mohr, 
Mathiesen and Evrard 1999; Ettori and Fabian 1999; Arnaud and Evrard 1999). A possible 
explanation of the ICM distributions is that the ICM has been heated non-gravitationally. In fact, 
Ponman, Cannon and Navarro (1999) indicated that the entropy of the ICM in the central regions 
of low-temperature or less massive clusters is higher than can be explained by gravitational collapse 
alone, although it is not understood what heats the ICM (e.g. supernovae or AGN) and where the 
ICM is heated, that is, outside or inside clusters. Heating other than the gravity of a cluster makes 
the ICM distribution flatter and different from the dark matter distribution. Thus, we expect that 
the position where a galaxy suffers from ram-pressure stripping depends on whether the ICM of 
the cluster (or the gas accreted by the cluster later on) has been heated non-gravitationally or not. 
In particular, we expect that the position where ram-pressure stripping occurs is more sensitive to 
the non-gravitational heating in the past. This is because a cluster breaks into progenitors or less 
massive clusters as time goes backwards, and because the heat required to explain the observations 
should have more influence on the ICM distributions of the less massive clusters (Cavaliere, Menci 
and Tozzi 1998; Balogh, Babul and Patton 1999; Loewenstein 2000; Fujita and Takahara 2000). 
Therefore, ram-pressure stripping in the progenitors may tell us when the ICM of present-day 
clusters was heated non-gravitationally; it will be a clue to the heating sources. 

Our paper is organized as follows. In §2 we describe our models of the dark matter distribution 
and the ICM distribution within clusters, the ram-pressure stripping of a radially infalling galaxy, 
and the evolution of cluster progenitors. In §3 we give the results of our calculations. We compare 
them with observations in §4. Conclusions are given in §5. 
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2. Models 

2.1. Distributions of Gravitational Matter and ICM 

The virial radius of a cluster with virial mass M v ; r is defined as 

3 M v - N 1/3 



rvir -UvrA c (.)p crit (z)J ' (1) 

where p cr it{z) is the critical density of the universe and A c (z) is the ratio of the average density of 
the cluster to the critical density at redshift z. The former is given by 

, x Pcrit,0^0(l +zf 

= W) , (2) 

where p C rit,o is the critical density at z = 0, and Q(z) is the cosmological density parameter. The 
latter is given by 

A c (z) = 18vr 2 (3) 

for the Einstein-de Sitter Universe and 

A c (z) = 18 7T 2 + 82x - 39x 2 , (4) 

for the flat universe with non-zero cosmological constant (Bryan and Norman 1998). In equation 
(4), the parameter x is given by x = Q(z) — 1. 

We assume that a cluster is spherically symmetric and the density distribution of gravitational 
matter is 

PmO) = Pmv(r/r vil )~ a , (5) 

where p mv and a are constants, and r is the distance from the cluster center. The normalization, 
Pmv, is given by 

3 — a 

Pmv = 3 A e p crit . (6) 

We choose a = 2.4, because the slope is consistent with observations (Horner, Mushotzky and 
Scharf 1999). Moreover, the results of numerical simulations show that the mass distribution in 
the outer region of clusters is approximately given by equation (5) with a ~ 2.4 (Navarro, Frenk 
and White 1996, 1997). 

We consider two ICM mass distributions. One follows equation (5) except for the normalization 
and the core structure; 

[l + (r/r c ) 2 ]-/ 2 

PiCM(r) = PiCM,vi r[1 + (rvir/rc)2] _ a/2 • (7) 
The ICM mass within the virial radius of a cluster is 

/•TVir 

Micm = / 4nr 2 picM(r)dr . (8) 
Jo 
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The normalization /3icM,vir is determined by the relation /b = M\qm/M v - u . where /b is the gas or 
baryon fraction of the universe. This distribution corresponds to the case where the ICM is in pres- 
sure equilibrium with the gravity of the cluster and is not heated by anything other than the gravity. 
We introduce the core structure to avoid the divergence of gas density at r = and use r c /r Y - u = 0.1. 
We call this distribution the 'non-heated ICM distribution'. We use /b = 0.25(/i/0.5)~ 3 / 2 , where 
the present value of the Hubble constant is written as Hq = 100 h km s _1 Mpc -1 . The value of 
fb is the observed ICM mass fraction of high-temperature clusters (Mohr, Mathiesen and Evrard 
1999; Ettori and Fabian 1999; Arnaud and Evrard 1999), for which the effect of non-gravitational 
heating is expected to be small. 

However, as mentioned in §1, observations show that ICM is additionally heated non-gravitationally 
at least for nearby clusters. Thus, we also model the distribution of the heated ICM using the ob- 
served parameters of nearby clusters as follows. In this paper, we assume that the ICM had been 
heated before accreted by clusters. However, the distribution will qualitatively be the same even if 
the ICM is heated after accreted by clusters (see Loewenstein 2000). 

Following Balogh, Babul and Patton (1999), we define the adiabat Kq = P/pg, where P is the 
gas pressure, p g is its density, and 7 is a constant. If ICM had already been heated before accreted 
by a cluster, the entropy prevents the gas from collapsing into the cluster with dark matter. In this 
case, the ICM fraction of the cluster is given by 



/iCM = min 



( M vir \ (Ko_y 3/2 ( t(z) 



(9) 



where -K34 = 10 34 erg g 5 / 3 cm 2 (Balogh, Babul and Patton 1999). 

The virial temperature of a cluster is given by 

k B T v i T 1 GM vir 



(10) 



fimu 2 r vir 

where k-Q is the Boltzmann constant, p,{= 0.61) is the mean molecular weight, tojj is the hydrogen 
mass, and G is the gravitational constant. When the virial temperature of a cluster is much 
larger than that of the gas accreted by the cluster, a shock forms near the virial radius of the 
cluster (Takizawa and Mineshige 1998; Cavaliere, Menci and Tozzi 1998). The temperature of the 
postshock gas (T 2 ) is related to that of the preshock gas (T\) and is approximately given by 

T 2 = T vir + |ri (11) 

(Cavaliere, Menci and Tozzi 1998). Since the gas temperature does not change very much for 
r < r v ; r (Takizawa and Mineshige 1998; see also Irwin and Bregman 2000), the ICM temperature 
of the cluster is given by Ticm = ^2- Since we assume that the density profile of gravitational 
matter is given by equation (5) with a = 2.4, the density profile of ICM is given by 

PiCM(r) = PiCM,v ir[1 + (rvir/rc)2] _ 3/3/2 , (12) 
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where (3 = (2.4/3)T vir /Ti C M (see Bahcall and Lubin 1994). Observations show that T x ~ 0.5-1 keV 
although it depends on the distribution of the gravitational matter in a cluster (Cavaliere, Menci and 
Tozzi 1998; Fujita and Takahara 2000). We choose 371/2 = 0.8 keV, hereafter. The normalization 
/9iCM,vir is determined by the relation /icm = Mccm/M™-- 

When T v i r < (3/2)Ti, a shock does not form but the gas accreted by a cluster adiabatically 
falls into the cluster. The ICM profile for r < r v ; r is obtained by solving the equation of hydrostatic 
equilibrium, 

dP GM(r) 

~T = 9 PICM , (.13) 

dr r z 

where M(r) is the mass of the cluster within radius r. Generally, equation (13) does not have 
analytical solutions for the matter distribution we adopted (equation [5] with a = 2.4). Thus, we 
use the solution for the isothermal distribution (p m [r] = /W,iso[?'Avir]~ 2 ) as an approximation; 
assuming 7 = 5/3, it is given by 



PICm(0 = PlCM,vir 

where 



1 + 



3 , /r vir \ 
A ln {—) 



3/2 

(14) 



A = (15) 
87rGp mv ,i so < ir 

(Balogh, Babul and Patton 1999). The parameter picM.vir is determined by the relation /icm = 
Micu/M VlT . In §3, we use the profile (12) for T vir > 3T x /2 (= 0.8 keV) and the profile (14) 
for T v ; r < 3Ti/2. We refer to this ICM distribution as 'the heated ICM distribution'. From 
observations, we use the value of Kq = O.37K34, which is assumed to be independent of cluster 
mass (Balogh, Babul and Patton 1999). 



2.2. A Radially Infalling Galaxy 

We consider a radially infalling disk galaxy accreted by a cluster with dark matter. The initial 
velocity of the model galaxy, Vi, is given at r = r v j r , and it is 

vj_ = GM vil _ GM vir (ig) 

2 ?Vir r ta 

where r ta is the turnaround radius of the cluster. Assuming that r ta = 2r v ; r on the basis of the 
virial theorem, the initial velocity is 

v * = y—- (17) 

The virial radius r v ; r is given by equation (1). 

The velocity of the model galaxy is obtained by solving the equation of motion; 

d 2 r GM(r) 

^2 " J5- • ( 18 ) 
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As the velocity of the galaxy increases, the ram-pressure from ICM also increases. The condition 
of ram-pressure stripping is 



PlCMVrel > 2"7tG£*£hI 



2.1 x 10- n dyncm- 2 ( ^-^r) 



10 kpc/ \8 x 10 20 mH cm -2 



(19) 



where £* is the gravitational surface mass density, Shi is the surface density of the HI gas, v ro t 
is the rotation velocity, and R is the characteristic radius of the galaxy (Fujita and Nagashima 
1999; Gunn and Gott 1972). Abadi, Moore and Bower (1999) have numerically confirmed that 
this analytic relation provides a good approximation. We define the cluster radius at which the 
condition (19) is satisfied for the first time as the stripping radius, r st . Since we assume that the 
ICM is nearly in pressure equilibrium for r < r v ; r , the relative velocity v rc \ is equivalent to the 
velocity of the galaxy relative to the cluster for r < r v ; r . 

Since the mass distribution of a cluster within the virial radius does not change rapidly, the 
time that a galaxy takes from r v ; r to r st is approximately given by 

At = / - , (20) 

Jr °* Jvf - 2A(f)(r) 

Where 4nGo r a 

A0(r) = ( , P ™ ™ , (r 2 ~ Q - r^) . (21) 
(2 — a){6 — a) 

Thus, when the ISM of the galaxy that was located at r = r v i r (i') at t = t' is stripped at r = r st at 
t = t' + At, the virial radius of the cluster becomes larger than r v i r (0- 



2.3. The Growth of Clusters 



We also investigate ram-pressure stripping in progenitors of a cluster. In this subsection, we 
construct a model of the growth of clusters. 

The conditional probability that a particle which resides in a object ('halo') of mass M2 at 
time t-2 is contained in a smaller halo of mass M\ ~ M\ + dM± at time t± {t\ < t^) is 



P 1 {M 1 ,t 1 \M 2 ,t 2 )dM 1 



1 



3c2 



2^ K - a 2 2 )3/2 



dMi 



cxp 



(S, 



cl 



s c2 ) 2 



2{a\ 



(To 



dMi , 



(22) 



where S c i is the critical density threshold for a spherical perturbation to collapse by the time tj, 
and <7j[= a (Mi)] is the rms density fluctuation smoothed over a region of mass Mj for i = 1 and 2 
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(Bond et al. 1991; Bower 1991; Lacey and Cole 1993). In this paper, we use an approximative 
formula of <5 c j (Nakamura and Suto 1997) and an fitting formula of a(M) for the CDM fluctuation 
spectrum (Kitayama 1998). 

We define the typical mass of halos at t that become part of a larger halo of mass Mq at later 
time io(> i) as 



M(t\Mo, to) = J ~ ~ ,„ |,„ ' ; ... , (23) 



J^ in P 1 (M,t\M ,t )dM 

where M m i n is the lower cutoff mass. We choose M m i n = 10 8 M Q , which corresponds to the mass 
of dwarf galaxies. 

In §3.2, we investigate a cluster progenitor whose virial mass is given by 

M viT (t\M , t ) = M(t\M ,t ) . (24) 

In that subsection, we will often represent M v j r (t|Mo, to) with M v [ T (t) or M v i r unless it is misun- 
derstood, and we will often abbreviate a cluster progenitor to a 'cluster'. 



3. Results 

The parameters for a model galaxy are v rot = 220 km s -1 , R = 10 kpc, and Ehi = 8 x 
10 20 mH cm~ 2 . Although the values are those for our Galaxy (Spitzer 1978; Binney and Tremaine 
1987), the condition of ram-pressure stripping (relation [19]) is not much different even in the 
case of smaller galaxies. Taking M33 for instance, the relation (19) turns to be picM^rei > 2-2 x 
10~ n dyn cm~ 2 , because the values are v ro t = 105 kms" 1 , R = 5 kpc, and Shi = 18 x 10 20 mH cm -2 
(Newton 1980; Garnett et al. 1997; Sofue et al. 1999). 

As cosmological models, we consider a standard cold dark mater model (SCDM model) and 
a cold dark matter model with non-zero cosmological constant (ACDM model). The cosmological 
parameters are h = 0.5, Qq = 1, Ao = A/(3Hq) = 0, and as = 0.63 for the SCDM model, and 
h = 0.75, 9,o = 0.25, A = 0.75, and a 8 = 1.0 for the ACDM model. As will be seen, the results do 
not much depend on the cosmological models. 



3.1. The Relation Between Mass and Stripping Radius 

We first discuss the results of the non-heated ICM models. Figure 1 shows the relations 
between the stripping radius, r st , and virial mass, M vir , of clusters at several redshifts. For a given 
redshift, r st is an increasing function of M V - 1T . One of the reasons is the mass-dependence of r v i r 
(Figure 2). That is, massive clusters are large. To see the effect of ram-pressure stripping, we 
illustrate r st /r v i r in Figure 3. For a given redshift, r st /r v i r is larger for more massive clusters. Since 
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not all galaxies have the pure radial orbits in real clusters, we suppose that in the real clusters with 
larger r st /r v ; r , more galaxies are affected by ram-pressure stripping. In that sense, ram-pressure 
stripping is more effective in more massive clusters. The mass dependence can be explicitly shown 
as follows. For a given redshift, the virial density of clusters, p v ir(= A c p cr i t ), does not depend on 
mass (equations [2] -[4]). Thus, the virial radius of a cluster is given by r vir oc (M vir / p vir y/ 3 oc mV 3 . 
Therefore the velocity of a galaxy infalling into the cluster has the relation, vf oc 

-^vir/'Vir °C -^yir " 

Since the typical ICM density of a cluster is proportional to the virial density, the ram-pressure of 
the galaxy follows the relation, picMvf oc M^/ 3 . 

The virial radius, r v i r , of a cluster with a given mass is smaller at higher redshift (Figure 2), 
because the virial density of the cluster p v ; r is increases as a function of redshift (equations [l]-[4]). 
On the other hand, r s t of a cluster with a given mass is larger at higher redshift in the non-heated 
ICM models (Figure 1). This is because the typical ICM density of clusters increases as the virial 
density increases. Since ram-pressure is proportional to ICM density, it is more effective at higher 
redshift (Figure 3). This can be clarified easily in the SCDM and non-heated ICM model as follows. 
The ICM density of clusters follows picm oc p vir oc (1 + z) 3 . For a fixed mass, the virial radius 
is represented by r v - lT oc (M v ; r / p v i T ) 1 ^ 3 oc (1 + z)^ 1 . Thus, the velocity of a galaxy is given by 
vf oc M v i r /r v ; r oc 1 + z, and the ram-pressure is given by picMvf oc (1 + z) 4 . This explains large 
r s t/r v ir of high-redshift clusters (Figure 3). 

Next, we show the results of the heated ICM models. In comparison with the non- heated ICM 
models, r st and r st /r v i r for less massive clusters at a given redshift in the heated ICM models are 
small (Figures 1 and 3). This is because of the small ICM density of the model clusters. As the 
mass of clusters decreases, (3 in equation (12) decreases, then the ICM distribution changes from 
equation (12) to (14), then the ICM fraction changes from /icm = fb to /icm < fb (equation [9]). 
In Figures 1 and 3, the small jump at M vir = 1.3 x 1O 14 M and z = in the ACDM model is due to 
the shift from equation (12) to (14). Figure 3 shows that in the heated ICM model, ram-pressure 
stripping is ineffective for galaxy clusters (or groups) with M v ; r < 5 x 10 13 M at z ~ and for 
those with M v ; r < 10 13 M at z ~ 1 — 2. However, for rich clusters with M v j r > 10 14 M observed 
at z ~ 0.5 — 1, the ram-pressure stripping should be more effective than the clusters with the same 
mass at z ~ regardless of the non-gravitational heating. 



3.2. The History of Ram-Pressure Stripping in Clusters 

Using the model constructed in §2.3, we investigate the history of ram-pressure stripping of 
galaxies in clusters. Figure 4 shows the evolutions of M v ; r (t|Mo, to) for Mq = 2 x 10 15 M and 
2 x 10 14 M , that is, the evolutions of clusters with typical mass. The redshifts corresponding to 
to are zq = and 0.5. We show the evolutions of r v - lv and r st in Figures 5 and 6, respectively. 
At z = zq, the stripping radius r st in the model of the heated ICM is not much different from 
that in the model of the non-heated ICM. In fact, the non-gravitational heating does not affect 
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ram-pressure stripping for clusters with mass of M v ; r > 10 14 M Q at z ~ (Figure 3). For clusters 
in that mass range, even if the gas had been heated before accreted by clusters, the clusters have 
gathered a large amount of gas and their virial temperatures have become large enough until z = zq. 
In our model of the heated ICM, these respectively mean that /icm = fb (equation [9]) and the 
ICM distribution is given by equation (12). At z = zq, the values of Picm are not much different 
between the heated ICM distribution and the non-heated ICM distribution at r ~ r st , although the 
slope of the former distribution is smaller. 

As z increases, r st in the SCDM model decreases faster than that in the ACDM model (Fig- 
ure 6). This is because M v - ir and r v ; r decrease faster in the former model (Figures 4 and 5). 
Moreover, r st decreases faster in the model of the heated ICM than that in the model of the non- 
heated ICM. In order to see the effect of ICM heating rather than that of the decrease of cluster 
size, we show the evolutions of r st /r v \ r in Figure 7. As can be seen, r st /r v ; r decreases rapidly at 
high redshifts in the models of the heated ICM, while it does not change significantly in the models 
of the non-heated ICM. The changes of the slope in the former models are due to the shift in the 
ICM distribution from equation (12) to (14) and the shift in the ICM fraction from /icm = fb to 
/icm < fb (equation [9]). Taking the heated ICM model of M = 2 x 10 15 M , z = 0, and SCDM 
as an example (Figure 6a and 7a), the shift in the ICM distribution occurs at z = 0.8 and the shift 
in /icm occurs at z = 1.0. The rapid decrease of r st /r v ; r in the model of the heated ICM is chiefly 
attributed to the decrease of /icm- On the other hand, the almost constant r st /r v i r in the models of 
the non-heated ICM is explained by the fact that although the mass of cluster progenitors and the 
velocity of galaxies in them decrease with z (Figure 4), the average mass density of progenitors, Pvir? 
and thus the average ICM density of progenitors, fbPvh, increase (equations [2]-[4]). We take the 
model of M = 2 x 10 15 M , z = 0, and SCDM as an example again. Figures 4 and 5 respectively 
show that M v ; r oc (1 + z)~ 6 and r v - 1T oc (1 + z)~ 3 . Thus, the velocity of the model galaxy has the 
redshift dependence such as vf oc M v ; r /r v i r oc (1 + z)~ 3 . On the other hand, picm oc p v - ir oc (1 + z) 3 . 
Therefore, the ram-pressure is almost independent of redshift or picuvf oc constant. 

In summary, Figure 7 shows that if the ICM (or the gas accreted by a cluster later on) is 
heated non-gravitationally at z » 1, ram-pressure stripping does not occur in cluster progenitors 
at z > 1 — 2. On the other hand, if the ICM had not been heated non-gravitationally until z ~ 0, 
ram-pressure stripping occurs even at z ~ 3. The difference can be explained by the small mass 
of cluster progenitors. At z ~ 1 — 2, the masses of progenitors are M v - ir < 10 13 M (Figure 4). 
Figure 3 shows that the non-gravitational heating significantly reduces the effect of ram-pressure 
at this mass range at z ~ 1 — 2. 

4. Discussion 

In the previous sections, we have modeled the ram-pressure stripping in clusters and their 
progenitors. In this section, we compare the results with several observations. Since the direct 
observation of ram-pressure stripping of galaxies, such as the observation of HI distribution in 
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galaxies, is difficult except for nearby clusters at present, we discuss the morphology and color of 
galaxies. 

In §3.1, we show that the ram-pressure stripping in clusters should be more effective in more 
massive clusters for a given redshift (Figure 3). Edge & Stewart (1991) investigate observational 
data of nearby clusters (z < 0.1) and found that the spiral fraction decreases and the SO frac- 
tion increases with the X-ray temperature and luminosity (see also Bahcall 1977; McHardy 1978). 
The clusters with high X-ray temperature and luminosity are generally massive ones (e.g. Horner, 
Mushotzky and Scharf 1999). Thus, the relations confirmed by Edge & Stewart (1991) are consis- 
tent with our prediction, if the ram-pressure stripping converts spiral galaxies into SO galaxies. 

Recently, clusters at z ~ 1 are observed in detail (Postman, Lubin and Oke 1998; Kajisawa 
et al. 2000). Lubin et al. (1998) observed a cluster CL 1604 + 4304 at z = 0.90 and found that the 
fraction of SO galaxies is large. The ratio of SO galaxies to elliptical galaxies in the cluster is 1.7+0.9, 
which is larger than that in intermediate redshift clusters (z ~ 0.5). On the contrary, the fraction 
of spiral galaxies is only 24^24%- Lubin et al. (1998) estimated that the virial mass of the cluster is 
M v ; r = 7.8^2'! x 10 14 /i _1 M Q , which seems to be exceptionally massive at that redshift. Our model 
predicts that at z > 1, most of the galaxies infalling into clusters of the mass should be subject to 
ram-pressure stripping (Figure 3). Thus, the high SO fraction may be due to the transformation 
of the field spiral galaxies by ram-pressure stripping. If this is the case, the transformation by 
stripping must be rapid. This is because the infall rate of field galaxies, most of which are blue 
spiral galaxies, increases with z (Kauffmann 1995); the ram-pressure stripping must convert rather 
part of the blue spiral galaxies into SO galaxies in a short time. On the other hand, most of the 
galaxies in a poor cluster CL 0023 + 0423 at z = 0.84 are normal spiral galaxies (Lubin et al. 1998). 
The population suggests that the ram-pressure stripping is not effective in the cluster. Since this 
cluster consists of two small components of M v - lT = l.Oto'4 x 10 13 /i _1 M Q and 2.6^Qg/?. _1 x 10 14 M , 
it is qualitatively consistent with our prediction (Figure 3). More observations of galaxy groups 
with M vir < 3 x 10 13 M at z ~ 1 may give us information about the non-gravitational heating 
of ICM for z > l(see Figure 3). However, for the quantitative comparison between the theory 
and the observations, we may need to use a so-called semi-analytic model of cluster formation or a 
numerical simulation including the effects of the galaxy infall rate, the variation of galaxy orbits, 
and the ram-pressure stripping. 

Tanaka et al. (2000) found that UV-excess red galaxies are abundant in a rich cluster near the 
quasar B2 1335 + 28 at z = 1.086; the fraction of such galaxies is relatively small at z ~ 0.2 (Smail et 
al. 1998). Tanaka et al. (2000) estimated the Abell richness of the cluster and found that it is class 
one. By comparing the Abell catalogue (Abell 1958) with X-ray catalogues (e.g. Jones and Forman 
1999; Mohr, Mathiesen and Evrard 1999; Ettori and Fabian 1999; Arnaud and Evrard 1999), it 
is shown that the ICM temperature of the class one clusters is Ticm ~ 3 — 4 keV. If Ticm = T v i r , 
the temperature corresponds to M V - 1T ~ 5 — 10 x 10 14 M Q in our models. Thus, our model predicts 
that most of the galaxies infalling into the cluster should be subject to ram-pressure stripping 
(Figure 3). Since the color of the UV-excess red galaxies can be explained by the superposition of 
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very weak star formation activity on old stellar population, the galaxies may be the ones suffering 
from ram-pressure stripping and the very weak star formation may be ember although alternative 
interpretations may be possible. If this is true, the abundance of the galaxies may also indicate 
that ram-pressure stripping happens very often at z ~ 1, although more sophisticated models are 
required for quantitative arguments as is the case of CL 1604 + 4304 and CL 0023 + 0423. In the 
future, it will be useful to investigate the morphology of the UV-excess red galaxies in high-redshift 
clusters in order to know whether the galaxies are the ones subject to ram-pressure stripping; if the 
galaxies have red disks, they are probably the ones (Fujita and Nagashima 1999). It is to be noted 
that in a rich cluster MS 1054-03 at z = 0.83, most of the spiral galaxies are red (van Dokkum et 
al. 2000). 

At intermediate redshift 0.2 < z < 0.5, observations show that the fraction of blue galaxies 
in rich clusters, /&, is larger than that at z ~ (Butcher and Oemler 1978, 1984). On the other 
hand, our model predicts that the fraction of galaxies affected by ram-pressure stripping increases 
with z. It seems that these contradict to each other, because the ram-pressure stripping suppresses 
star formation in galaxies (Fujita and Nagashima 1999). However, as is mentioned above, the 
infall rate of field spiral galaxies is higher at the intermediate redshift in comparison with z ~ 0. 
Thus, the increase of fb that occurs in galaxy clusters at the redshift may suggest that the effect of 
the increase of the infall rate overwhelms the effect of the ram-pressure stripping in snapshots of 
galaxy population, and that the increase of fb is not inconsistent with our prediction; we predict 
that the for a given mass, ram-pressure stripping at the intermediate redshift is not as effective as 
that at z ~ 1 (Figure 3). Moreover, the average mass of the clusters observed at the intermediate 
redshift is expect to be smaller than that at higher redshift (z ~ 1), because the clusters at the 
intermediate redshift are nearer to us and easier to be observed. On the other hand, observations 
also suggest that the star formation of the infalling galaxies are ultimately truncated in clusters 
at the intermediate redshift, although there is a contradiction over the way the star formation 
is suppressed (Poggianti et al. 1999; Balogh et al. 1999). Considering that most of the clusters 
observed at the intermediate redshift are fairly rich, ram-pressure stripping should be fairly effective 
in the clusters, although it is not as effective as that at z ~ 1 (Figure 3). Thus, our model appears 
to be consistent with the truncation of star formation. Moreover, it is qualitatively consistent with 
the observed transformation of spiral galaxies into SO galaxies in rich clusters at the intermediate 
redshift (see §1). 

Except for extremely massive clusters (M vir > 3 x 10 15 M ), Fi gure 1 shows that r st < 1 Mpc 
regardless of redshift z and the ICM heating. For nearby clusters at z ~ 0, observations show that 
the fraction of SO galaxies increases at r < 2 Mpc (Whitmore, Gilmore and Jones 1993; Dressier 
et al. 1997), which appears to be larger than the prediction by our ram-pressure model if the 
transformation into SO galaxies is due to ram-pressure stripping. Moreover, for clusters observed at 
0-2 ^ z < 0.5, the fraction of blue galaxies also decreases at r < 2 Mpc (Abraham et al. 1996; van 
Dokkum et al. 1998; Morris et al. 1998; Kodama and Bower 2000). If ram-pressure stripping is 
the main mechanism of the transformation of blue spiral galaxies into SO galaxies, the observations 
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suggest that some galaxies at r > 1 Mpc have already been affected by the stripping when they 
were in the progenitors of the clusters before accreted by the main cluster progenitors. In fact, 
Figure 7 shows that the ram-pressure stripping in progenitors is effective at least for z < 3 if the 
ICM is not heated non-gravitationally and it is effective for z < 1 — 2 even if the ICM is heated non- 
gravitationally. Moreover, even if galaxies were inside the stripping radius of the main progenitor 
at some earlier time, some of them may have been scattered to large apocenter orbits during the 
merger process of cluster progenitors (Fujita et al. 1999; Balogh, Navarro and Morris 2000). 

Analyzing observational data of galaxies in rich clusters, Kodama and Bower (2000) investi- 
gated the star formation history of all the galaxies in the central regions of the clusters. They found 
that star the formation rate per galaxy mass declines more rapidly than in the field environment 
at z < 1; it suggests the truncation of star formation in most of the galaxies. This may imply that 
ram-pressure stripping has been effective in the clusters or in their progenitors at least for z < 1, 
because the star formation rate of galaxies should decline after the ISM, from which stars are born, 
is stripped. Since Figure 7 shows that the ram-pressure stripping has been effective for z < 1 for 
rich clusters regardless of the non-gravitational heating, the results of Kodama and Bower (2000) 
are consistent with our predictions. Unfortunately, because of large uncertainty, their results can- 
not constrain the star formation rate of the galaxies in the cluster progenitors for z > 1. Thus, we 
cannot discuss the effect of the non-gravitational heating (Figure 7). 



5. Conclusions 

We have studied ram-pressure stripping of galaxies in clusters and their progenitors. In partic- 
ular, we pay attention to its dependence on redshift and the mass of clusters. As a model galaxy, 
we consider a radially infalling disk galaxy; the initial position and velocity are given by a spherical 
collapse model of structure formation. Since X-ray observations show that the ICM of nearby clus- 
ters is heated non-gravitationally, we also investigate the effect of the heating on the ram-pressure 
stripping. Our main findings are the following: 

1. For a given redshift, ram-pressure stripping of galaxies is more effective in more massive 
clusters. This is because the velocity of the radially infalling galaxy increases with the virial mass 
of the cluster. If ram-pressure stripping transforms spiral galaxies into SO galaxies, our model 
is consistent with the observed relation between galaxy populations and cluster luminosities (or 
temperatures) . 

2. For a given mass of clusters, ram-pressure stripping of galaxies in the clusters is more 
effective at higher redshift. This is because the density of the intracluster medium increases with 
the redshift. In particular, at z > 1, most of the galaxies radially infalling into the centers of rich 
clusters are affected by ram-pressure stripping. The relatively high fraction of SO galaxies and the 
abundance of UV-excess galaxies in rich clusters at z ~ 1 may be due to the ram-pressure stripping. 

3. The non-gravitational heating reduces the effect of ram-pressure stripping for clusters with 
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-Mvir < 5 x 10 13 M at z ~ and for those with M vir < 10 13 M at z ~ 1. However, for clusters 
with M v i r > 10 14 Mq, it does not have an influence on the effect of ram-pressure stripping. 

4. If the ICM (or the gas accreted by a cluster later on) is heated non-gravitationally at 
z » 1, ram-pressure stripping does not occur in cluster progenitors at z > 1 — 2, because the heat 
makes the ICM fraction of the cluster progenitors small. On the other hand, if the ICM had not 
been heated non-gravitationally until z ~ 0, ram-pressure stripping occurs even at z ~ 3. 

I am grateful to T. Yamada, M. Nagashima, I. Tanaka, T. Kodama, T. Tsuchiya, and D. A. 
Dale for useful discussions and comments. Comments from an anonymous referee led to significant 
improvements in the quality of this paper. 

REFERENCES 

Abell, G. O. 1958, ApJS, 3, 211 

Abadi, M. C, Moore, B. and Bower, R. G. 1999, MNRAS, 308, 947 

Abraham, R. G. and 8 colleagues 1996, ApJ, 471, 694 

Arnaud, M. and Evrard, A. E. 1999, MNRAS, 305, 631 

Bahcall, N. A. 1977, ApJ, 218, L93 

Bahcall, N. A. and Lubin, L. M. 1994, ApJ, 426, 513 

Balogh, M. L., Babul, A. and Patton, D. R. 1999, MNRAS, 307, 463 

Balogh, M. L., Morris, S. L., Yee, H. K. C, Carlberg, R. G. and Ellingson, E. 1999, ApJ, 527, 54 
Balogh, M. L., Navarro, J. F. and Morris, S. L. 2000, ApJ, 540, 113 
Bekki, K. 1998, ApJ, 502, L133 

Binney, J. and Tremaine, S. 1987, Galactic Dynamics (Princeton:; Princeton Univ. Press) 

Bond, J. R., Cole, S., Efstathiou, G. and Kaiser, N. 1991, ApJ, 379, 440 

Bower, R. G. 1991, MNRAS, 248, 332 

Bryan, G. L. and Norman, M. L. 1998, ApJ, 495, 80 

Burstein, D. 1979, ApJ, 234, 435 

Burstein, D. 1979, ApJ, 234, 829 

Butcher, H. and Oemler, A. 1978, ApJ, 219, 18 



-15- 



Butcher, H. and Oemler, A. 1984, ApJ, 285, 426 

Byrd, G. and Valtonen, M. 1990, ApJ, 350, 89 

Cavaliere, A., Menci, N. and Tozzi, P. 1998, ApJ, 501, 493 

Couch, W. J., Barger, A. J., Smail, I., Ellis, R. S. and Sharpies, R. M. 1998, ApJ, 497, 188 
Dale, D. A. & Uson, J. M. 2000, AJ, 120, 552 

Dressier, A., Oemler, A. J., Butcher, H. R. and Gunn, J. E. 1994, ApJ, 430, 107 
Dressier, A. and 8 colleagues 1997, ApJ, 490, 577 
Edge, A. C. k Stewart, G. C. 1991, MNRAS, 252, 428 

Ellingson, E., Lin, H., Yee, H. K. C. and Carlberg, R. G. 2000, ApJ, in press 

Ettori, S. and Fabian, A. C. 1999, MNRAS, 305, 834 

Fujita, Y. 1998, ApJ, 509, 587 

Fujita, Y. and Nagashima, M. 1999, ApJ, 516, 619 

Fujita, Y. and Takahara, F. 2000, ApJ, 536, 523 

Fujita, Y., Takizawa, M., Nagashima, M. and Enoki, M. 1999, PASJ, 51, LI 
Gaetz, T. J., Salpeter, E. E. and Shaviv, G. 1987, ApJ, 316, 530 

Garnett, D. R., Shields, G. A., Skillman, E. D., Sagan, S. P. and Dufour, R. J. 1997, ApJ, 489, 63 
Giovanelli, R. & Haynes, M. P. 1985, ApJ, 292, 404 
Gunn, J. E. and Gott, J. R. I. 1972, ApJ, 176, 1 

Horner, D. J., Mushotzky, R. F. and Scharf, C. A. 1999, ApJ, 520, 78 

Irwin, J. A. and Bregman, J. N. 2000, ApJ, 538, 543 

Jones, C. and Forman, W. 1999, ApJ, 511, 65 

Kajisawa, M. and 19 colleagues 2000, PASJ, 52, 61 

Kauffmann, G. 1995, MNRAS, 274, 153 

Kitayama, T. 1998, PhD thesis, Univ. Tokyo. 

Kodama, T. and Bower, R. G. 2000, MNRAS, in press (astro-ph/0005397) 
Lacey, C. and Cole, S. 1993, MNRAS, 262, 627 



-16- 



Loewenstein, M. 2000, ApJ, 532, 17 
Lubin, L. M. 1996, AJ, 112, 23 

Lubin, L. M., Postman, M., Oke, J. B., Ratnatunga, K. U., Gunn, J. E., Hoessel, J. G. & Schneider, 
D. P. 1998, AJ, 116, 584 

McHardy, I. 1978, MNRAS, 184, 783 

Mohr, J. J., Mathiesen, B. and Evrard, A. E. 1999, ApJ, 517, 627 

Moore, B., Katz, N., Lake, G., Dressier, A. and Oemler, A. 1996, Nature, 379, 613 

Mori, M. and Burkert, A. 2000, ApJ, 538, 559 

Morris, S. L., Hutchings, J. B., Carlberg, R. G., Yee, H. K. C, Ellingson, E., Balogh, M. L., 
Abraham, R. G. and Smecker-Hane, T. A. 1998, ApJ, 507, 84 

Nakamura, T. T. and Suto, Y. 1997, Prog. Theor. Phys., 97, 49 

Navarro, J. F., Frenk, C S. and White, S. D. M. 1996, ApJ, 462, 563 

Navarro, J. F., Frenk, C. S. and White, S. D. M. 1997, ApJ, 490, 493 

Newton, K. 1980, MNRAS, 190, 689 

Poggianti, B. M., Smail, I., Dressier, A., Couch, W. J., Barger, A. J., Butcher, H., Ellis, R. S. and 
Oemler, A. J. 1999, ApJ, 518, 576 

Ponman, T. J., Cannon, D. B. and Navarro, J. F. 1999, Nature, 397, 135 

Portnoy, D., Pistinner, S. and Shaviv, G. 1993, ApJS, 86, 95 

Postman, M., Lubin, L. M. and Oke, J. B. 1998, AJ, 116, 560 

Rakos, K. D. and Schombert, J. M. 1995, ApJ, 439, 47 

Smail, I., Edge, A. C, Ellis, R. S. and Blandford, R. D. 1998, MNRAS, 293, 124 

Sofue, Y., Tutui, Y., Honma, M., Tomita, A., Takamiya, T., Koda, J. and Takeda, Y. 1999, ApJ, 
523, 136 

Spitzer, L. 1978, Phisical Processes in the Interstellar Medium (A Wiley- Interscience Publication:; 
New York) 

Quilis, V., Moore, B. and Bower, R. 2000, Science, 288, 1617 

Takeda, H., Nulsen, P. E. J. and Fabian, A. C. 1984, MNRAS, 208, 261 

Takizawa, M. and Mineshige, S. 1998, ApJ, 499, 82 



-17- 



Tanaka, I., Yamada, T., Aragon-Salamanca, A., Kodama, T., Miyaji, T., Ohta, K. & Arimoto, N. 
2000, ApJ, 528, 123 

van Dokkum, P. G., Franx, M., Kelson, D. D., Illingworth, G. D., Fisher, D. and Fabricant, D. 
1998, ApJ, 500, 714 

van Dokkum, P. G., Franx, M., Fabricant, D., Illingworth, G. D. & Kelson, D. D. 2000, ApJ, 541, 
95 

Whitmore, B. C., Gilmore, D. M. and Jones, C. 1993, ApJ, 407, 489 



This preprint was prepared with the AAS IATjrjX macros v5.0. 



-18- 




Fig. 1. — The relations between the virial mass and the stripping radius of clusters at z = 0, 1, 
and 2. (a) SCDM, (b) ACDM. The solid and dotted lines indicate the heated and the non-heated 
ICM models, respectively. 




Fig. 2. — The relations between the virial mass and the virial radius of clusters at z = 0, 1, and 2. 
(a) SCDM, (b) ACDM. 
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Fig. 3. — The relations between the virial mass and the ratio of the stripping radius to the virial 
radius at z = 0, 1, and 2. (a) SCDM, (b) ACDM. The solid and dotted lines indicate the heated 
and the non-heated ICM models, respectively. 




Fig. 4. — The evolutions of the virial mass of clusters, (a) zq = 0, (b) zq = 0.5. Four models are 
plotted in each figure; SCDM and M = 2 x 10 15 M (thin solid), SCDM and M = 2 x 10 14 M 
(thin dotted), ACDM and M = 2 x 10 15 M© (thick solid), and ACDM and M = 2 x 10 14 M© 
(thick dotted). 
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Fig. 5. — The evolutions of the virial radius of progenitors, (a) zq = 0, (b) z$ = 0.5. Four models 
are plotted in each figure; SCDM and M = 2 x 10 15 M (thin solid), SCDM and M = 2 x 10 14 M 
(thin dotted), ACDM and M = 2 x 10 15 M© (thick solid), and ACDM and M = 2 x 10 14 M Q 
(thick dotted). 
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Fig. 6. — The evolutions of the stripping radius, (a) SCDM and zq = 0, (b) SCDM and zq = 
0.5. (c) ACDM and zq = (d) ACDM and zq = 0.5. Four models are plotted in each figure; 
M = 2 x 10 15 M and the non-heated ICM distribution (thin solid), M Q = 2 x 10 14 M and the 
non- heated ICM distribution (thin dotted), M = 2 x 10 15 M and the heated ICM distribution 
(thick solid), and M = 2 x 10 14 M and the heated ICM distribution (thick dotted). 
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Fig. 7. — The evolutions of the ratio of the stripping radius to the virial radius, (a) SCDM and 
z = 0, (b) SCDM and z = 0.5. (c) ACDM and z = (d) ACDM and z = 0.5. Four models 
are plotted in each figure; Mq = 2 x 10 15 M and the non- heated ICM distribution (thin solid), 
M = 2 x 10 14 M and the non-heated ICM distribution (thin dotted), M = 2 x 10 15 M and the 
heated ICM distribution (thick solid), and M = 2 x 10 14 M and the heated ICM distribution 
(thick dotted). 



